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Abstract 

Two-particle correlations in invariant mass are studied separately for like-sign and unlikc-sign 
charged particles produced in deep inelastic positron-proton scattering in a new kinematical 
domain. The data were taken with the HI detector at the HERA storage ring in 1994, in 
which 27.5 GeV positrons collided with 820 GeV protons at a centre of mass energy y/s = 300 
GcV. 

The observed enhancement of the like-sign correlations at low invariant masses is related 
to the dimensions of the hadronic source. The data are compared to different QCD models 
where the hadronization is performed with the string-fragmentation model. Results are 
presented for the first time separately for diffractive and non-diffractive scattering, in domains 
of four-momentum transfer, Bjorken-cc, and hadronic center of mass energy, and in intervals 
of charged particle multiplicity. The observed source radii do not differ strongly from those 
measured in lower energy lepton-nucleon inelastic scattering, and e^e^ annihilation. 
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1 Introduction 



Intensity interferometry was first used in radio astronomy by Hanbury-Brown and Twiss (HBT- 
effect H]) to measure the spatial extent of stars in 1953. In particle physics the HBT-effect was 
rediscovered in 1959 by G. and S. Goldhaber, Lee and Pais in antiproton-proton annihilation Q. 
The GGLP effect, as it was called then, in its original form, described the tendency of pairs of 
like-sign pions to have smaller opening angles than pairs of unlike-sign pions. This effect was 
shown to arise from the symmetrisation of the wave functions of multi-boson states, hence is 
also known as Bose-Einstein correlation, and now normally refers to the enhanced probability 
for identical boson pairs to have similar momenta. The shape of the enhancement in relative 
momentum space was shown to be related to the spatial dimensions of the pion source. Bose- 
Einstein correlations (BEG) have since been studied in high-energy physics in nearly every type 
of reaction, e.g. in e^e~ annihilation, in lepton nucleon (^A^) scattering, in vrA^, KN , NN, and 
NN reactions, as well as in relativistic heavy-ion collisions with the aim of understanding the 
space-time evolution of the boson source. Many reviews of this subject exist @-|6|. In recent 
years it has been emphasized that Bose-Einstein correlations are but one of the many facets of 
multi-particle correlation studies (see e.g. ||^^ for a review). Independently of the details of the 
interpretation, the experimental study of BEG in a broad variety of reactions should eventually 
lead to better understanding of the space-time geometry of multihadron production 



1.1 Bose-Einstein correlations 



Two-particle Bose-Einstein correlations are usually described in terms of normalized two-particle 
densities 

P(1)P(2) ' ^ ^ 

where P(l,2) and P{i) {i = 1,2) are, respectively, the two-particle and single-particle inclusive 
densities. Bose-Einstein correlations show up as an enhancement in the production rate of iden- 
tical pions with similar four-momenta, and reflect both geometrical and dynamical properties of 
the particle radiating source. If the particles are assumed to be emitted from a set of indepen- 
dently radiating sources distributed in space-time with a distribution p{Cj, and the particles are 
identical bosons, imposing Bose symmetry on the production amplitude leads to a correlation 
function 

i?(r) = i?o(l + A I /5(T) |2) , (2) 

where p is the Fourier transform of p{Cj] Rq is a normalisation constant and T is related to the 
invariant mass M of a pion pair with four-momenta pi and p2 by 

= -{pi - P2f = - 4ml . (3) 



The parameter A in Eq. (|2|) is usually assumed to reflect the degree of incoherence of the 
source, i.e. A = 1 for complete incoherence and A = for complete coherence. A static sphere 
of emitters with Gaussian density of the form 



p(O = p(0)exp(-|^) , (4) 
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where the parameter r corresponds to the size of the production volume, leads to a correlation 
function 

fi(T) = i?o(l + Aexp(-r2r2)) , (5) 

usually referred to as Goldhaber parametrisation 

However, the simple geometrical interpretation of the interference pattern, based on an 
optical analogy is invalid when emitters move relativistically with respect to each other, leading 
to strong correlations between the space-time and momentum-energy coordinates of emitted 



particles |TC, 11]. Correlations of this type arise due to the nature of inside-outside cascade 



dynamics |l^ as in colour-string fragmentation |13]. In the interpretation of BEG by Andersson 



and Hofmann |14] and Bowler 1 11,15] in the string model, the length scale measured by BEG is 
therefore not related to the size of the total pion emitting source, but to the space-time separation 
between production points for which the momentum distributions still overlap. This distance is, 
in turn, related to the string tension. The model predicts an approximately exponential shape 
of the correlation function 

ii(T) = i?o(l + A exp(-rr)) , (6) 
where r is expected to be independent of the total interaction energy. 

In recent years, considerable experimental effort has been devoted to the study of fluctu- 
ation phenomena in multiparticle production processes These studies, complemented by 
theoretical analyses of the (approximately) self-similar character of perturbative QGD cascades, 
concentrated on the search for "intermittency" : the occurrence of large fluctuations in particle 
density arising from scale-invariant dynamics. This work has provided firm indications that 
intermittency is strongly connected with Bose-Einstein correlations Scale invariance implies 
that multiparticle correlation functions exhibit power law behaviour over a considerable range 
of the relevant relative distance measure (such as or M^) in phase space. As such, BEG from 
a static source do not exhibit power law behaviour. However, a power law is obtained if the 
size of the particle source fluctuates event-by-event, and/or, if the source itself is a self-similar 
(fractal-like) object extending over a large volume ]jl6| . Power law behaviour of BEG has been 
observed in several experiments p7|-p^. In these studies, the ratio R is parametrised using the 
form 

R{M) = A + B(^-^y . (7) 

1.2 Outline 

In this paper we analyse BEG in deep inelastic positron proton scattering (DIS) at HERA. The 
present study differs from previous studies of BEG in iN scattering p0|-[2^ in several aspects. 
Since the data were taken at a tenfold higher centre of mass energy, they also cover a larger range 
in the kinematical quantities Bjorken— x and momentum transfer ]|2^. The correlations in 
the data can therefore be analysed separately as a function of x, Q^, and hadronic centre of 
mass energy W. In addition diffractive and non-diffractive events can be distinguished. 

At HERA about 10% of the deep inelastic scattering events exhibit a large region of rapidity 



in which there is no significant hadronic activity ]]24,25]. These events are usually attributed to 
diffractive processes and are well described by Monte Garlo simulations in which diffraction is 
modeled as deep inelastic positron pomeron scattering where the pomeron has partonic structure. 
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To the best of our knowledge, BEC have never been studied in diffractive processes. Although it 
is expected that the fragmentation mechanism in these events is similar to that in non-diffractive 
DIS, it nevertheless remains of interest to analyse BEC separately for that class of events. 

Aside from comparing diffractive and non-diffractive samples, a specific goal of our analysis is 
to confront the correlation data with the different interpretations referred to above. In addition 
we use the Goldhaber parametrisation as a tool for the comparison with previous experiments. 

This paper is organized as follows. After a brief presentation of the HI detector, the event 
and track selections are described in Sect. ^. The Monte Carlo models used for comparison are 
listed in Sect. ^ The details of the analysis procedure and systematic uncertainties are given in 
Sect. ^. The results and conclusions are presented in Sect. |6| and Sect. |^ respectively. 



2 The HI detector 

The data were taken with the HI detector at the HERA ep - collider in 1994, with 27.5 GeV 
positrons colliding head-on with 820 GeV protons at a centre of mass energy ^/s = 300 GeV. 
A detailed description of the HI detector can be found elsewhere Here we briefly describe 
the components relevant for this study. 

The present analysis is based on the central tracking system, closest to the interaction point, 
and covering polar angles between 20° and 160°, where 9 is defined with respect to the proton 
beam direction {+z axis). Measurements of the charge and momentum of charged particles are 
provided by two large cylindrical drift chambers (central jet chambers, CJCl and CJC2). The 
chambers have wires strung parallel to the beam axis (z— direction) with the drift cells inclined 
with respect to the radial direction, yielding up to 56 space points. The measured space point 
resolution in the drift coordinate {r — (f) plane) is 170 ^m and can, by comparing signals read 
out at both wire ends, achieve a resolution of one percent of the wire length in z. Two sets of 
thin cylindrical z— drift chambers complement the measurement of charged track momenta and 
are located at the inner radii of the jet chambers. The z— chambers provide track elements with 
typically 300 /im resolution in z and 8° in (p. This requires a drift direction parallel to, and 
sense wires perpendicular to, the beam axis. The central tracking system is complemented by a 
forward tracking system and a backward proportional chamber (BPC) covering the polar angle 
ranges 7° < 9 < 25° and 155° < 9 < 175°, respectively. The BPC combined with a measurement 
of the interaction vertex position, is used to measure the polar angle of the scattered positron 
to better than 1 mrad. 

Surrounding the tracking system is a liquid argon calorimeter consisting of electromagnetic 
and hadronic sections covering 4° < < 153° over the full azimuth. Outside these detectors a 
superconducting coil provides a uniform magnetic field of 1.15 T parallel to the beam line. To 
measure the energy of the scattered positron, a backward electromagnetic calorimeter (BEMC) 
covers the angular range 155.5° < 9^ < 174.5°. Behind the BEMC is the time-of-flight system 
(ToF), consisting of two planes of scintillators, with a time resolution of 1 ns. It is used to reject 
events caused by the interaction of the proton beam with material before the detector. Such 
particles arrive at the ToF system before particles from the nominal interaction point. In the 
forward direction, surrounding the beam-pipe, is the plug, a copper /silicon calorimeter, covering 
the angular range 0.7° < 9 < 3.5°. Further forward are the forward muon drift chambers on 
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either side of a toroidal magnet around the beam-pipe. Their primary purpose is to detect low 
angle muons, but they are also used to detect particles coming from secondary interactions with 
collimators designed to protect the detector from synchrotron radiation. 

The most relevant performance parameters of the HI detector for this analysis are the in- 
variant mass resolution for charged particle pairs and the resolution of reconstructed event 
kinematics (for the latter see ref. |2^). The two-particle invariant mass resolution can be de- 
duced from the widths of the signal at M = 0.5 GeV and of the J/^ peak at M = 3.1 GeV 
to be au = 8.5 MeV and au = 76 MeV respectively. This agrees with Monte Carlo 
studies which also predict ar = 3.5 MeV at r=0.1 GeV. 



3 Event and track selection 

The present study is based on a sample of low DIS events. The data sample comprises an 
integrated luminosity of 1.26 pb~^. The scattered positron can be detected in the BEMC for 
in the range from 1 GeV^ to 100 GeV^. The event selection used in this analysis follows closely 
that described in our structure function analysis [ 23 1 . In brief, the trigger used to collect low 



data required a localised energy deposit in the BEMC of at least 4 GeV and no veto from the 
ToF system. This has been shown to be more than 99% efficient for the selection of DIS events 
in which the scattered positron entering the BEMC has an energy above 10 GeV. 



3.1 Event selection 

All events are required to have a vertex, reconstructed from tracks in the central tracker, located 
within lb 30 cm of the mean vertex z-position. This cut is required to reject beam-induced 
background and to permit a reliable determination of the kinematic variables. To exclude events 
with large QED radiative effects and to ensure substantial hadronic energy flow in the detector 
we require the invariant mass squared, W"^, of the hadronic system determined both with the 



Jacquet-Blondel method |29| and from the scattered positron to be larger than 4400 GeV^, as well 
as the scattered positron to have an energy E^, > 12 GeV and a polar angle 157° < 6e < 172.5°. 
The energy constraint ensures that the remaining photoproduction background is less than 
1%. Finally, at least two selected tracks are required in each event irrespective of charge. 
This selection accepts events for in the range 6 < < 100 GeV^ and for x in the range 
10"'* <x< 10~^. Both X and are calculated from the positron nergy and angle. 

The data are divided into diffractive and non-diffractive subsamples. Non-diffr active events 
are selected by requiring that the energy deposited in the liquid argon calorimeter within the 
polar angular range 4.4° < 9 < 15° be greater than 0.5 GeV. Diffractive events are selected by 
demanding that the energy deposited in the plug is less than 3 GeV, that there are fewer than 
three pairs of hits in the forward muon drift chambers, and that the most forward significant 
energy deposit (above 400 MeV) detected in the liquid argon calorimeter has a pseudo-rapidity 
of 3.0 or less. A further cut is made using the variable Xp, defined as 

Q2 + Ml 



X 



(8) 



where is the mass of the hadronic final state separated from the proton remnant system by 
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the largest rapidity gap 1 24 1 . In order to restrict the data to the kinematic region where pomeron 
exchange dominates and to ensure good acceptance of the dissociating pomeron-virtual boson 
system it is required that < 0.05. The quantity Xp can be interpreted as the fraction of the 
proton momentum carried by the pomeron. 

After all cuts, the non-diffractive sample contains 48 000 events and the diffractive sample 
2 500 events. The same selections are applied to the various Monte Carlo generated event samples 
used in the analysis. Further details are given in references [30,31|. 



3.2 Track selection 

Charged tracks are required to be linked to the primary event vertex and to show at least 10 
hits within the jet chambers. Requiring a transverse momentum pt > 0.15 GeV rejects tracks 
curving strongly within the CJC. A polar angle 9 satisfying 22° < 9 < 150° restricts the analysis 
to particles which can be measured in both jet drift chambers in order to improve the quality 
of the tracks. 

BEC studies are sensitive to double counting effects, i.e. the splitting of a single track into 
two pieces, which mimic two like-charge tracks with a very low value of T. In this experiment 
we find that tracks are split mainly between the inner and outer jet chamber. Thus only tracks 
which start within the inner CJC are considered. Monte Carlo studies have shown that this 
selection has a 95% track finding efficiency, and the remaining rate of split tracks is 0.3 %. 
This was verified by a visual scan of data. In addition no spurious tracks due to noise hits 
or reconstructed mirror hits of real tracks were found. To remove sensitivity to the remaining 
rate of split tracks all fits to the data are performed in the region T > 0.018 GeV. To further 
assess the quality of the simulation program, a comparison of track parameters between data 
and reconstructed Monte Carlo was made. Pairs with low invariant masses were selected both 
in data and in Monte Carlo and technical quantities such as radial track length, which is a 
good measure of track quality, were compared. Good agreement was found down to the lowest 
invariant mass included in this analysis. 

All charged particles are assumed to be pions in this analysis. The effect of this assumption is 



discussed in Sect. 5.2. The influence of small differences in the reconstruction efficiency between 



positive and negative tracks due to the geometry of the jet chamber and a small torsion and 



consequent misalignment of the two jet chambers |32] were taken into account by separately 
analysing pairs of both charges. The effect on the results is described in Sect. ^.2| and found to 
be small. 



4 Monte Carlo models 

Non-diffractive neutral current DIS events are modelled using the Monte Carlo generator LEPTO 
|3^,|3^ which uses exact QCD matrix elements to first order to describe boson-gluon fusion, 
and QCD-Compton processes. Soft parton radiation is described by either the colour dipole 
model (CDM) |^5| as implemented by ARIADNE 4.03 |3^, in which the partons are radiated 
independently from a chain of dipoles, or the parton shower approach |37|, where soft parton 



radiation is generated in the "leading log" approximation. The combination of matrix elements 
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with matched parton showers is referred to as the MEPS model pj], whilst combining matrix 
elements with CDM is referred to as MEAR and uses LEPTO 6.1 



To study the effects of QED radiation the generator DJANGO 6.0 |38| is used which is an 
interface between HERACLES 4.4 and LEPTO. It allows the inclusion of 0{a) processes 
including real photon emission from the positron as well as virtual electro-weak corrections. 



The Monte Carlo generator RAPGAP 2.1 |4C] is used to model diffractive DIS. The model 



uses parton density functions taken from leading order QCD fits to the measured diffractive 
structure function F^^^^ As with LEPTO, exact matrix elements are used to calculate 

first order QCD radiation, with additional radiation incorporated using either CDM, or parton 
showers. In addition RAPGAP includes virtual QED effects to first order and real photon 
emission. The model employing RAPGAP with CDM is referred to as RAPA, and RAPGAP 
with parton showers is referred to as RAPP (Fig. ^). 

All Monte Carlo generators used in this analysis employ the Lund string hadronisation model 
as implemented in JETSET 7.4 [42|. A method for the simulation of Bose-Einstein correlations 



is available as an option in the JETSET Monte Carlo |43]. The algorithm does not rigorously 
model BEC by symmetrising production amplitudes, but treats it as a final state interaction 
by reshuffling hadron momenta. As an option, the model assumes a Gaussian shape of the 
Bose-Einstein enhancement^ according to the Goldhaber prescription of Eq. (|^ . This was used 
to generate a sample of non-diffractive DIS events referred to as MEAR(BEC). Unless stated 
otherwise, all Monte Carlo models refer to those without BEC simulated. 



The HI detector simulation program, HISIM, is based on the CERN package GEANT |44| 
and was used to simulate the detector response in detail and to correct the data for geometrical 
acceptance, kinematical cuts, resolution and particle interactions with the detector material. 
The simulated events were processed through the same reconstruction chain as real data. 



5 Analysis procedure 

5.1 Construction of the correlation function 

In order to measure the correlation function of Eq. (||) we normalize the two-particle like-sign 
inclusive density denoted by P2{T) to a reference sample p^^^ {T) which ideally contains no 
Bose-Einstein correlations and form the ratio 

R(T) = ^f-y ^ith ^'2^-^^^^^ and j p{{T)dT^{n^) . (9) 

Here A'" is the total number of events, n^^{T) the number of (++) and ( ) pairs in the sample, 

and <v}> is the mean number of like-sign pairs per event. 

The choice of the reference sample is not trivial and has been a source of bias and systematic 
errors in all BEC measurements to date. Ideally it should satisfy the following conditions: 

^The parameters MSTJ{51) and MSTJ{52) in the generator code were changed from their default values 
(see [ p2|p3[ ) to 2 and 7, respectively. This choice leads to the inclusion of BEC for multi-pion and multi-kaon states 
with a Gaussian parametrisation, for which we have chosen PARJ{92) = A = 1 and PARJ{93) = (hc/r) — 0.38 
GeV, for r=0.53 fm. 
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absence of BEC, presence of correlations due to the topology and the global properties of the 
events, and absence of dynamical correlations due to resonances not present in like-sign pairs. 
In this analysis we use either the two-particle unlike-sign inclusive distribution 

p-/(r) = p«(r)^l^^^^P^ with |pf(T)dr^(n"), (lo) 

where n^^(T) is the number of (H — ) pairs in the sample, and (n") is the mean number of 
unlike-sign pairs per event, or we create uncorrelated pairs by mixing tracks from different 
events, denoted by 

pvef ^j,^ = p^(T) ^ normalized as j pi® pi{T) dT = (n') . (11) 

A mixed event is constructed by taking all combinations of each track in one event with all 
tracks in another event. In order to reduce statistical errors each event is mixed with 20 others. 

We then define three ratios of densities as follows: 

j:^/m P2 . . j^lu P2 (12) 

Plf^ Pi ' Pi® Pi ' P2 ' 



The same techniques were commonly used in the past in many experiments fl^,!!,!!,!!,!!-!!, 
but for each experiment special treatment is necessary. 

At HERA, events are observed in a wide range of virtual-photon proton (7*^) centre of 
mass energies, and hence the data sample comprises quite different topologies depending on the 
position of an event in the x — Q^ plane. Events with large on average have a current (struck 
quark) system which is boosted further forward in the laboratory system than events at lower 
Q^. Event topologies as seen in the transverse (r — (/)) plane have a hadronic system which is 
balanced against the scattered lepton. Mixing tracks from events with non-collinear hadronic 
systems, would yield track pairs with artificially large momentum differences. Thus, in order to 
overcome this bias, all events are rotated in the r — <j) plane such that the positron has = 0. 
A similar effect is expected to hold in the 9 — z plane. However, this effect can be reduced by 
mixing events with similar and W values. Mixed events are, therefore, required to differ in 

by not more than 20 GeV^, and in W by less than 20 GeV. The difference in event charged 
multiplicity was restricted to be less than three. 

The unlike-sign, like-sign and event-mixed densities constructed from our data are shown 
for the non-diffractive sample in Fig. ||a and for the Monte Carlo simulations in Fig. |l]b. Since 
the standard Monte Carlo simulation does not include BEC, the like-sign pair distribution 
constitutes, in principle, the ideal reference distribution for a Monte Carlo with BEC included 
(shown as stars in Fig. P^). Neither the mixed-event nor the unlike-sign reference distribution 
from data is consistent with p2 obtained from the Monte Carlo model without BEC. 

Mixing tracks from different events apparently shifts the distribution towards larger T com- 
pared to the invariant mass distribution of the same event, as we can deduce by comparing 
Pi ®> Pi{T) with P2{T) , both shown in Fig. We therefore conclude that despite the pre- 
cautions taken in constructing it, the mixed-pair sample still differs from the ideal reference 
distribution in the topological, momentum and charge conservation constraints. 

The unlike-pair distribution p2 contains dynamical correlations from Kg and resonance de- 
cays, most notably from the p^, 77, r]', and uo. A sharp peak at T ~ 0.4 GeV can be seen in 
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P2 (Fig. |T]a) which arises from the Kg and the broad maximum at T ~ 0.7 GeV stems from 
the decay. The contributions from r/, ry', and w decay do not appear as sharp peaks and are 
more troublesome for the analysis than those from Kg and p^, because they contribute to the 
low mass region, where BEG appear. To illustrate this we have used the Monte Garlo model to 
display the contributions from resonance decays in Fig. The K^ (not shown) and p^ decays 
cluster in a sufficiently narrow mass region, and are removed by excluding the regions 0.37 — 0.43 
GeV (K^) and 0.65 - 0.85 GeV (/9°) in T from the analysis. 

The ratio of the data to the MEAR(BEG) model shown in Fig. ||c decreases smoothly with 
T for both p2 and P2- Distinct dips or peaks are not observed in either the low mass, the Kg 
or the p^ region. A smooth T dependence is also observed for the ratio of the density for pion 
pairs to the density for all charged particle pairs. This ratio was obtained within the MEAR 
Monte Carlo model and is shown in Fig. ^3. 

The shortcomings of the reference distributions constructed from the data can be corrected, 
to first order, if both the correlated and uncorrelated distributions are divided by the corre- 
sponding Monte Carlo distributions. Thus in order to discriminate BEG from other dynamical 
correlations, we consider a double ratio by dividing the correlation function of Eq. (P) obtained 
from the data, by that obtained from reconstructed Monte Carlo events which do not contain 
BEG: 

= j^McJ^ . (13) 

This procedure should also correct for the detector acceptance, analysis cuts and the lack 
of particle identification provided that the Monte Garlo simulation can represent the observed 



distributions sufficiently well, as we will show in Sect. 6.1. The double ratio is then fitted with a 
modified version of Eq. (|5|) or Eq. (^) in order to extract the parameters r and A for the Gaussian 
and the exponential parametrisations: 

RR{T)= i?o(l + ar)(l + Aexp(-r2r2)) , (14) 
RR{T)= i?o(l + ar)(l + Aexp(-rr)) . (15) 

Here Rq is a global normalization factor, and the additional parameter a allows for remaining 
long-range correlations. 



5.2 Systematic errors and corrections 

To verify the stability of our results with respect to the assumptions made concerning the refer- 
ence distributions, and to estimate the systematic errors arising from experimental resolutions 
and the treatment of the data, we have performed a number of checks described in the following. 
In order to quantitatively assess the influence of an assumption or a particular analysis step, 
the deviations of the parameter values relative to that of the reference fit are determined. The 
reference fit uses Eq. ([l^) for the Gaussian parametrisation, and Eq. (|l5|) and Eq. (^ for the 
exponential and power law fits, respectively. The MEAR Monte Garlo is used to calculate the 
double ratios, and the fit to the non-diffractive data is performed in the range from 0.018 to 2 
GeV in T (0.018 to 1.0 GeV in T for power law fits) and to the diffractive data in the range 
from 0.036 to 2 GeV. For the unlike-sign reference distribution, the regions of the p^ resonance 
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and the are excluded from the fit (see Sect. The observed deviations from the reference 
fit are then considered in the evaluation of the systematic error. 

We first discuss those checks which use Monte Carlo events only, and where we study the 
influence of reference sample choice and detector effects. 

Reference sample: The systematic effects associated with the choice of the reference sam- 
ple can be studied "in Monte Carlo only" experiments using MEAR(BEC) treated as 
"data". MEAR without BEC is used for the denominator of the double ratios. This 
can be done at both the reconstructed level, i.e. after full detector simulation and re- 
construction, and at the generator level. In addition, the like-sign density from MEAR 
without BEC may be used ideally in the denominator for RiT) in Eq. (^, with p2 in the 
numerator obtained from the MEAR(BEC) Monte Carlo model. In this case the single 
ratio /92(MEAR(BEC))/p2(MEAR) is fitted with the standard Gaussian parametrisation 
(Eq. (p!^)). The results are given in Table |l] and show that a fit using the ideal reference 
distribution is able to reproduce the input value for r (0.53 fm) within the statistical errors. 
The results of fits to the double ratios RR^"^ and iii?'" are also given in Table |l|. Inspection 



of these results reveals, not unexpectedly (see Sect. |5T| above), that neither the use of the 
unlike-sign nor the event-mixed reference distribution can reproduce Monte Carlo input 
values, not even at the generator level. The shift towards larger radii, between 0.18 fm 
for unlike-sign pairs and 0.08 fm for event-mixing is significant, and could, in principle, be 
used to correct the experimental results. We have not applied such a correction to allow 
comparison with other experiments (Sect. |6.3| ) which have observed this shift, but have 
not corrected for it. 

The value extracted for A (0.38) differs from the input value (A = 1), since resonance 
decay products and charged particles other then pions are not distinguished when con- 
structing the correlations. It is with respect to this value, that we must evaluate further 
differences introduced by the choice of the reference sample. Table |^ shows that these are 
sizeable. For this reason, the values of A found by different experiments, and in different 
reactions, differ significantly and it seems hazardous to draw definitive conclusions from 
such comparisons. 

Detector Effects: Finite detector resolution and acceptance as well as imperfect pattern recog- 
nition in the reconstruction step can influence the BEC analysis. As Table |^ shows the 
difference in the extracted parameters between generator and detector level for the radius 
parameter is small, though clearly noticeable for A. This comparison is repeated for all sets 
of low non-diffractive data discussed below separately and included as a contribution 
to the systematic error. 

The following contributions to the systematic error are evaluated from the data directly. 

Models: To assess the model dependency of the results we utilize both the MEAR and MEPS 
Monte Carlo models for the construction of the acceptance corrected double ratios. 

Background parametrisation: As alternatives to the linear parametrisation of long-range 
correlations used in our reference fits (Eq. ( p!^ ) and Eq. (|l5|)) we used the following forms 

RR{T) = Ro{l + aT + eT^){l + Xexp{-r^T^)) , (16) 
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Monte 
Carlo 
Level 


Reference samples 


like-sign P2iT) 
r(fm) A 


unlike-sign RR^^ 
r(fm) A 


event-mixed RR^"^ 
r(fm) A 


GEN 


0.51 ± 0.02 0.38 ± 0.02 






REC 
GEN 




0.67 ± 0.04 0.48 ± 0.04 
0.71 ± 0.04 0.57 ± 0.04 


0.61 ± 0.04 0.28 ± 0.02 
0.61 ± 0.03 0.35 ± 0.02 



Table 1: Gaussian fit results (Eq. (|T|)) using MEAR(BEC) Monte Carlo (see Sect. |) treated as 
"data". The input parameters for MEAR(BEC) were r = 0.53 fm, and A = 1. MEAR without 
BEC is used as the denominator for construction of the double and single ratios. The entries 
"GEN" and "REC" refer to using the Monte Carlo before and after the detector simulation and 
reconstruction step. Results from fits to the single ratio are given in the column labeled 
Results from fits to the double ratios RR^^ and RR^^ are given in the last two columns. 



RR{T)= i?o(l + aT + er2)(l + A exp(-rr)) . (17) 

Assuming that the uncorrelated background is constant over the full range in T, as in 
Eq. (|5|) and Eq. (|6|), leads to worse values of in most cases, and hence is not considered. 
We also repeated all Gaussian and exponential fits with an upper limit in T of 1.2 GeV 
instead of 2 GeV. For the power law fit (Eq. (|^)) the upper limit was increased to 1.2 GeV. 

Charge sign: Differences in the detector response to positively and negatively charged particles 
could cause systematic shifts of the data. Therefore we analysed positive and negative pairs 
separately to assess the influence of such effects. For the diffractive data we have assumed 
the same error as for the non-diffractive sample, since further splitting of an already 
statistically small sample would lead to statistical errors far exceeding the systematic 
effect we were studying. 

Data treatment: Technical changes such as bin width, track selection criteria, and the influ- 
ence of the resonances (including the and mass region for the unlike-sign reference 
sample) were investigated and found to have no influence outside statistical fluctuations. 
In addition the restriction T > 0.018 GeV was raised to 0.036 GeV and led to negligible 
changes in all fit parameters. 

Event-mixing: The analysis was repeated using a simpler method of event-mixing whereby the 
data sample was subdivided into four W classes and each event was mixed with one other 
event within the same class. In this analysis a less restrictive W"^ selection (as determined 
via the Jacquet-Blondel method) was used, in addition to an alternative track selection. 



The two analyses are in good agreement with each other [30|. Furthermore, the analysis of 
the diffractive data was repeated with more stringent restrictions on the mixing of events. 
Specifically, the requirements were tightened to ISW < 10 GeV, and AQ^ < 10 GeV^. 
This was also found to have no influence on the results pi]]. 

Final-state interactions: We have not applied corrections for electromagnetic |50, 51] and 
strong ]52] vrvr interactions in the final state to the data, or for the purity of the pion 



sample, as e.g. done by the DELPHI-collaboration jl7]. We calculated the Coulomb 
corrections and found that they affect only the region T < 0.05 GeV, changing A and r by 
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6X = 0.03 and 6r = 0.02 fm, respectively. A purity correction in the sense, that one tries 
to correct the observed distribution eUminating decay products from long-Uved resonances 
and charged particles other than pions, can be made by rescaling the data using the ratio 
of pions, which are either directly produced or decay products of short lived resonances, 
to that of all charged particles given by the Monte Carlo models. Using this procedure we 



have confirmed the observations of the DELPHI-collaboration [ 17 1 , that only changes in 
the radius parameter of about la (statistical) are observed, whilst A becomes consistent 
with one (see 0,0)- 

QED radiation: QED effects were investigated using the Monte Carlo model DJANGO p8| 



and found to be negligible for the ratios used in the analysis |31]. 



The most important contributions to the systematic errors for diffractive and non-diffractive 
data are given in Table ^ and Table ^, and are added in quadrature for the final errors quoted. 
For the exponential and power law parametrisations (Eq. ( p!5[ ) and Eq. (|^) the corresponding 
systematic uncertainties are listed in Table ^ and Table |5[ When using the preferred event-mixed 
technique, the systematic errors of r are sufficiently small to allow a meaningful comparison 
with other experiments and between different subsets of our data. However, since systematic 
uncertainties in the determination of A are as large as 60% in some cases, and A is also very 
sensitive to pion purity, we place no emphasis on our measurement of A, but include results in 
the following for completeness. 





unlike-si^ 


;n p!(t) 


event-mixed pi pi (T) 




(5r(fm) 


6X 


5r(fm) 


6X 


Detector effects 


-0.04 


-0.09 


±0.00 


-0.06 


Background parametrisation 


-0.01 


±0.01 


±0.02 


±0.00 


Models 


±0.02 


±0.19 


±0.00 


±0.06 


Track charge 


±0.01 


±0.01 


±0.02 


±0.02 


Sum 


±0.02 


±0.19 


±0.03 


±0.06 


(quadratic) 


-0.05 


-0.21 


-0.02 


-0.09 


Reference sample f 


±0.18 


±0.19 


±0.08 


-0.03 


Statistical error 


±0.04 


±0.03 


±0.03 


±0.02 



Table 2: Contributions to the systematic errors for the non-diffractive data sample and the 
Gaussian parametrisation (Eq. (|T^)). The last row shows the statistical error obtained from the 
reference fit for comparison. ^ Observed shift in the Monte Carlo extracted parameters compared 
to a like-sign reference distribution without BEC (see text and Table Q). 



6 Results 

6.1 Comparison of the data with Monte Carlo predictions 

In Figs. ^, Q and|5| we examine how well the data are described by the Monte Carlo simulations. 
From the ratio i?"™(T) shown in Fig. ^a for data, MEAR and MEPS Monte Carlo we can infer 
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unlike-si^ 




event-mixed 


pi (g> pi{T) 




(5r(fm) 


5X 


5r(fm) 


6X 


Background parametrisations 


±0.04 


±0.23 


±0.01 


±0.13 


Models 


±0.03 


±0.11 


±0.01 


±0.08 


Track charge 


±0.01 


±0.01 


±0.02 


±0.02 


Sum 


±0.05 


±0.26 


±0.02 


±0.15 


(quadratic) 


-0.05 


-0.11 


-0.03 


-0.08 


Statistical error 


±0.13 


±0.13 


±0.06 


±0.08 



Table 3: Contributions to the systematic errors for the diffractive data sample using the Gaussian 
parametrisation ([l^). The statistical error is that obtained from the reference fit. 





unlike-si 
6r{im) 


gn p^(T) 
6X 


event-mixed pi pi (T) 
(5r(fm) 5X 


Detector effects 


±0.02 


±0.15 


-0.04 


±0.11 


Background parametrisation 


-0.27 


±0.54 


±0.08 


-0.09 


Models 


±0.04 


±0.38 


±0.04 


±0.12 


Track charge 


±0.01 


±0.02 


±0.03 


±0.05 


Sum 


±0.05 


±0.68 


±0.09 


±0.17 


(quadratic) 


-0.27 


-0.38 


-0.06 


-0.16 


Statistical error 


±0.09 


±0.08 


±0.11 


±0.06 



Table 4: Contributions to the systematic errors for the non-diffractive data sample using an 



exponential parametrisation (15). The statistical error is that obtained from the reference fit. 





unlike-sign P2{T) 


event-mixed pi 5 






6(3 


SB 


SA 


6/3 


5B 


5A 


Detector effects 


-0.01 


-0.001 


-0.01 


-0.09 


±0.006 


-0.01 


Background parametrisation 


±0.29 


-0.002 


±0.02 


±0.12 


-0.005 


-0.01 


Models 


±0.30 


±0.003 


±0.03 


±0.08 


±0.003 


±0.00 


Track charge 


±0.15 


±0.003 


±0.00 


±0.04 


±0.001 


±0.00 


Sum 


±0.44 


±0.004 


±0.04 


±0.15 


±0.007 


±0.00 


(quadratic) 


-0.34 


-0.005 


-0.03 


-0.13 


-0.006 


-0.01 


Statistical error 


±0.20 


±0.002 


±0.01 


±0.15 


±0.006 


±0.01 



Table 5: Contributions to the systematic errors for the non-diffractive data sample using a power 
law parametrisation (|^). The statistical error is that obtained from the reference fit. 
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that the fragmentation model (JETSET) seems to overestimate the amount of production, 
as previously observed |17, 18|, |5|,^. However, as no such discrepancy is seen in Fig. ||c, we 
conclude that the inclusion of BEC in Monte Carlo models is able to satisfactorily describe the 
production region. The Monte Carlo also fails to predict correctly the T region (0.1 < T < 0.2 
GeV), where reflections of the tj, u and tj' resonances appear (Fig. ^). At the smallest values of 
T a sharp spike (truncated in the figures) can been seen in Fig. ^a and Fig. ^3. For like-sign pairs 
this is attributed to the residual effects of track double counting. Similarly, photons converting 
into e+e~ pairs contribute to the first bin in the unlike-sign pair distributions. Good agreement 
between data and Monte Carlo model is observed in the lowest bin (not visible in Figs. ^ and 

a 



In Figs, ^b, ^c, ^ and |^ the Bose-Einstein enhancement is clearly visible for like-sign pairs 
for T < 0.4 GeV in the ratios i?'™ and i?'". The slow rise of i?'™" towards threshold, predicted 
by the models, is caused by mixing tracks of events with a different topology as discussed in 



Sect. qJJ and in [|17|^. At larger T the Monte Carlo models give a good description of the 
data. The and signals are visible as distinct dips in Fig. 

In Fig. |5|a we compare the ratio i?"™ for non-diffractive data with the MEAR(BEC) Monte 
Carlo simulation. Apart from the very low T region, good agreement is seen. Figure ^b shows 
the ratio i?'™. Since the radius r used in the generation step lies close to the experimental 
result, this model provides a rather accurate representation of the data except for T < 0.2 
GeV where the data systematically exceed the prediction. This gives evidence that the Bose- 
Einstein enhancement rises somewhat faster towards threshold than expected from a Gaussian 
parametrisation of the correlation function. This was noticed in previous experiments (see 



e.g. |22, 4c,53,5l) 



6.2 Fits to the measured distributions 

6.2.1 BEC in diffractive and non-difFractive DIS 

Figures ^ and show the data, from which our final results are derived and which serve as 
reference plots for the discussion of the systematic errors in Sect. |5.2|. The results from the 
fits to these double ratios using the Goldhaber parametrisation (Eg. (|14[)) are superimposed on 
the data. Table |6| summarises the results, which were obtained in this analysis using the non- 
diffractive and diffractive data sets, for both the event-mixed and unlike-pair reference samples. 



The parameter a (typically between 0.01 to 0.02 GeV ) has been omitted (see Eqs. (14) - ([iq)). 

As seen from Table the BEC parameters derived from the event-mixed and from the unlike- 
sign reference samples differ substantially. This was also observed in many earlier experiments 
||l7,2C,^, H, 55, 56|, and subsequently ascribed @, |5^ to correlations due to the production 
of long-lived resonances and vrvr— interactions in the final state, as we have discussed above. 
This observation is independent of the choice of the Monte Carlo generator used to describe 
the scattering process. The difference appears already at the generator level, in the double and 
single ratios as is apparent from Table |^. The event-mixing method, however, leads to values 
of r which are in closer agreement with the input value used in the Monte Carlo generator with 
BEC. 

The data in Table ^ provide no evidence for a statistically significant difference between 
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Data set 


event-mixed pi ® Pi(T) 


r (fm) 


A 


xVndf 


non-diffractive 
diffr active 


0.54 ± 0.03 to^^ 

0.49 ± 0.06 j:[!:[]3 


0.32 ± 0.02 +^:^;^ 

0.46 ± 0.08 1[J;J| 


96/72 
18/23 


Data set 


unlike-sign P2{T) 


r (fm) 


A 


xVndf 


non-diffractive 
diffr active 


0.68 ± 0.04 to^l 
0.59 ± 0.13 


0.52 ± 0.03 
0.46 ± 0.13 t^'n 


77/56 
26/17 



Table 6: Summary of experimental results for the Gaussian fits (Eq. ([l^)) to the two data 
samples using both types of reference distributions. The first error is statistical, the second is 
systematic. 



the BEG parameters in predominantly non-diffractive and diffractive DIS interactions. The 
observation of similar BEG effects supports the idea that the basic fragmentation mechanism is 
similar for diffractive and non-diffractive events. 

When comparing with other experiments, or searching for dependencies on kinematic quan- 
tities such as rc, Q'^, W , or multiplicity, we concentrate in the following on the statistically more 
significant non-diffractive data set. 

6.2.2 Kinematical and multiplicity dependence of BEC 

In order to study BEC under different conditions, the non-diffractive data sample is divided into 
three classes for the three kinematic variables separately, as listed in Table 0. The double ratio 
correlation functions RR{T) are constructed with the mixed-event reference sample and then 
fitted to Eq. (0). The results are summarised in Table |7[ The parameters are found, within 
the statistical errors, to be independent of the kinematical region considered. 





r(fm) A 


r(fm) A 


r(fm) A 


X 


0.60±0.06 0.30±0.03 
(0.0001 <x < 0.0006) 


0.56±0.05 0.34±0.03 
(0.0006 <x< 0.0019) 


0.44±0.06 0.38±0.07 
(0.0019 <x< 0.01) 


Q' (GeV^) 


0.52±0.04 0.42±0.04 
(6 < Q2 < 12) 


0.63±0.08 0.25±0.04 
(12 < Q2 < 25) 


0.47±0.04 0.41±0.05 
(25 < Q2 < 100) 


W (GeV) 


0.52±0.07 0.26±0.05 

(65 < ^ < 120) 


0.48±0.03 0.42±0.04 
(120 <W < 180) 


0.68±0.08 0.34±0.04 
(180 <W < 240) 



Table 7: Parameters r and A extracted using Eq. ( [T^ ) for each subset of the non-diffractive data 
sample. Only statistical errors are given. 

We have also grouped the non-diffractive data into three subsets of observed charged particle 
multiplicity. Figure ^ shows the corresponding experimental distributions and fits. The classes 
and the corresponding values found for r and A are listed in Table ^. The value of the parameter 
r is observed to rise with increasing multiplicity. Indications for an opposite trend are noticeable 
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for X. Previously we have observed that the multiplicity is nearly independent of Q^, but depends 
weakly on W ||5^. The kinematical relation ~ xW"^, at small x, leads one then to expect a 
positive correlation of r with W and a negative correlation with x, if r increases with multiplicity. 
The data appear consistent with this expectation, though the effect is not significant in view of 
the statistical errors. 

For each bin of observed multiplicity we determine the corrected true mean multiplicity using 
the iterative matrix migration method of reference [^. The corrected mean charged particle 
densities are given in Table ^. The j*p centre of mass pseudo-rapidity interval, as determined 
by our track selection criteria, is Arj* = 3.2 (1.1 < r]* < 4.3). 

In order to assess the significance of the observed increase of r with multiplicity, we repeat 
the systematic error analysis described in Sect. |5.2| for each multiplicity bin. The influence of the 
shape of the background slope at large T is such that the systematic uncertainty (see Table ^) is 
twice as large as the statistical errors, and hence the rise consistent with zero. For comparison 
with recent e'^e~ data we have repeated the analysis with the unlike-sign reference sample, 
and also find a value for the slope statistically consistent with zero. 



Observed 
Multiplicity 


Corrected 
Multiplicity 


event-mixed pi (g) pi (T) 
r(fm) A 


unlike-sign P2{T) 
r(fm) A 


4 < n < 7 
7 < n < 12 
n > 12 


4.9 ± 1.1^ 
8.2 ± 1.6^ 
13.6 ± 2.4t 


0.42±0.05 0.37±0.05 
0.58±0.05 0.31±0.03 
0.81±0.12 0.42±0.07 


0.53±0.06 0.54±0.08 
0.77±0.07 0.54±0.06 
0.72±0.09 0.65±0.09 


Slope 


event-mixed pi (8) pi (T) 


unlike-sign P2{T) 


1 dr 

<r> dn 


0.085 ± 0.026 toml 


0.045 ± 0.023 t^^l 


1 dX 

<X> dn 


0.009 ± 0.018 toM3 


0.024 ± 0.026 t^^^l 



Table 8: Parameters r and A extracted using Eq. (|T^ for different multiplicity subsets of the 
non-diffractive data sample. The last two rows list the result for the relative slope of the two 
parameters with respect to charged multiplicity. The first column indicates the interval in 
observed multiplicity, the second column the corresponding corrected mean multiplicty. The 
first error quoted is statistical, and where given, the second is systematic. ^ These numbers are 
not errors, but represent the spread in true charged multiplicity in each subset. 



6.2.3 Alternative parametrisations of BEC 

Figure |9| shows the non-diffractive data of Fig. ^ together with the results from fits using the 
exponential, power law, and Gaussian parametrisations of Eq. (|l5|), Eq. (|^, and Eq. (p^), 
respectively. The parameters for the fits are listed in Table ^. The quality of the fit for the 
exponential is slightly better than for the Gaussian parametrisation (x^/ndf= 85/72 compared 
to 96/72) and confirms previous observations (see e.g. |Q) that the Bose-Einstein correlation 
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function is decreasing faster with T than a Gaussian. 

The power law form for the two-particle correlation function was originally suggested in the 
framework of a Mueller-Regge analysis of multiparticle inclusive spectra where the invariant 
mass occurs as a natural variable [Q. In recent studies of the intermittency effect, and its 
relation to BEG, power law behaviour for like-sign particle pairs was observed in hadron-hadron 
collisions and e^e~ annihilation at LEP [17| with values for the slope (3 close to 1.5. From 
Fig. I and Table | it is seen that a scale-invariant form of BEG provides a valid alternative 
parametrisation of the Bose-Einstein effect. In addition we find that fitting the power law over 
the range 0.018 < T < 2.0 GeV, and using the fit function 



RR{M) = A + €-M + B{ 



/3 



(18) 



which has an additional term linear in M with coefficient e as a free parameter, we are able 
to obtain good fits over the range 1 < T < 2 GeV, too. These result in a reduction of /? to 
0.91 ± 0.18 with xVndf=86/72 for RR^"", and /? = 1.49 ± 0.19 with xVndf=80/56 for RR^"". 



Power law 




B [GeV-^] 


A 


xVndf 


RR^"' 
RR^"" 


1.20 ± 0.15 to'il 
1.82 ± 0.20 t^f^ 


0.018 ± 0.006 1^:^^^ 
0.005 ± 0.002 1°:°°^ 


0.93 ± 0.01 IJ];^^ 
0.94 ± 0.01 tilt 


49/49 
52/33 


Exponential 


a [GeV-i] 


r [fm] 


A 


xVndf 


RR^'^ 
RR^"" 


0.08 ± 0.04 
0.13 ± 0.02 


0.68 ± 0.11 ta^^Q 
0.99 ± 0.09 tofj 


0.64 ± 0.06 
1.00 ± 0.08 to^^ 


85/72 
85/56 


Gaussian 


a [GeV-i] 


r [fm] 


A 


xVndf 


RR^"" 
RR^"" 


0.02 ± 0.01 
0.08 ± 0.02 


0.54 ± 0.03 
0.68 ± 0.04 taos 


0.32 ± 0.02 tl^^Q 
0.52 ± 0.03 t^Xi 


96/72 
77/56 



Table 9: Results from the fits to the double ratios RR^'^{T) and RR'''^{T) from the non-diffractive 
data with a Gaussian (Eq. (|l^)), exponential (Eq. (|l^)) and power law (Eq. (0)) parametrisation. 
Where given, the first error is statistical, the second systematical. 



6.3 Comparison with other experiments 



Table |10| summarises the Gaussian fit results from other experiments for comparison with the 
present analysis. For completeness also results are included, which used correction methods dif- 
ferent from ours, e.g. those where purity or Goulomb corrections were applied. These corrections 
mainly influence the parameter A, which may explain the larger differences observed here. The 
values of r extracted with a mixed-event reference sample are systematically lower than those 
obtained with an unlike-sign reference, consistent with our data. However, when allowance is 
made for systematic errors, which not all experiments have quoted, but which are likely to be of 
similar size to ours, the measured radii are relatively consistent with each other. The fact that 
the data from £A^— scattering are independent of the energy, and furthermore agree with those 
from e^e~ annihilation at various centre of mass energies (with the exception of the MARK II - 
data [^]), seems to lend support to the string model interpretation, where the source dimension 
is related to the string tension and is energy independent. 
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The multiplicity dependence has also been studied in a few experiments. In e^e^ annihilation 
at LEP a small rise of r with observed multiplicity has been reported by OPAL, while A 
was seen to decrease: 

1 dr 1 dX 

0.0036 ± 0.0006 , — ^ = - 0.0042 ± 0.0008 . 



< r > dn < X > dn 

The unlike-sign reference sample was used. Our values for this reference sample 

1 dr I dX 

0.045 ±0.023, — ^ - = 0.024 ± 0.026 , 



< r > dn < X > dn 

are consistent with these results, albeit within large errors. 

It may be argued that a more appropriate comparison is to examine the slopes with respect 
to mean particle density per interval of pseudorapidity rather then per interval of multiplicity, 
because the angular acceptance of the LEP experiments is different in the centre of mass system 
and the total multiplicity larger (10 < n < 40). We have scaled the e~^e~ results by the range 
in r]* accessible to OPAL (Arj* = 2 and | r/* |< 1) and also our data by Ar]* to obtain: 

1 UT" 1 (JiT' 

0.0072 ± 0.0012 (OPAL) , = 0.144 ± 0.074toMj (HI) . 



< r > d{dn / dr}*) ' < r > d{dn / drj*) 

The two data sets are still consistent. 

In proton-proton and proton-antiproton collisions, first indications for a multiplicity depen- 
dence were seen at low centre of mass energies at the ISR [^3|. Both the UAl-collaboration 
at centre of mass energies between 200 and 900 GeV (| rj* |< 3) as well as the E735-collabo- 
ration at 1.8 TeV (| rj* |< 3.25) report a strong increase of radius with particle density, 
dn/dr]* , in pp collisions. Figure ^ shows the results on r, together with our data converted into 
particle density, where all data were obtained using event-mixed reference samples. For UAl 
we have only plotted the statistically most significant data set at 630 GeV. Our results match 
these data at the lower end reasonably well. 



7 Conclusions 

We have presented a study of Bose-Einstein correlations (BEC) at the highest centre of mass 
energy for ep— collisions available so far. Even at this high energy (300 GeV) the apparent source 
dimensions, measured by the width of the BEC signal are in good agreement with data from iN 
scattering at lower energies (10 to 30 GeV) and with most data from e~^e~ reactions between 
^/s = 29 and ^/s = 91 GeV. BEC in diffractive DIS have been measured for the first time, and 
no difference between diffractive and non-diffractive data has been observed. No Q^, x oi W 
dependence is observed in the kinematical ranges accessible at HERA. These observations fit 
well within the picture inherent to string models, that hadronisation is a universal phenomenon 
and that the measured spatial dimensions reflect mainly the constant string tension. 

For the non-diffractive data, where statistics are highest, we have examined alternative 
parametrisations of the Bose-Einstein correlation function. Both an exponential form in T and 
a power law in invariant mass describe the Bose-Einstein enhancement at threshold well. The 
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Experiment 
Vi, W [GeV] 


unlike-sign P2{T) 
r{im) A 


event-mixed pi ® pi (T) 
r(fm) A 




0.83±0.03 0.31±0.02 
(0.82±0.03) (0.45±0.02) 

0.82±0.04 0.48±0.03 

0.955±0.012±0.015 0.672±0.013±0.024 

0.73±0.05±0.20 0.47±0.05±0.05 

0.82±0.06±0.04 0.35±0.03 
(0.80±0.06±0.04) (0.70±0.06±0.09) 
0.75±0.03±0.04 0.28±0.02±0.04 
(0.84±0.06±0.05) (0.50±0.03±0.04) 


0.47±0.03 0.24±0.02 
(0.42±0.04) (0.35±0.04) 
(0.49±0.01±0.05) (1.06±0.05±0.16) 
0.52±0.02 0.30±0.01 

0.58±0.06±0.02 0.39±0.04±0.03 

0.97±0.10±0.05 0.27±0.04±0.02 
(1.01±0.09±0.06) (0.45±0.03±0.04) 
0.65±0.04±0.05 0.50±0.04 
(0.65±0.04±0.05) (0.61±0.05±0.06) 


DELPHI m 

91 

DELPHI JitI 
ALEPH If 

Q1 

OPAL Hi 

91 

AMY 01 

58 

TASSO in 

34 

MARK lit [H 
29 

TPCt Hi 

29 


m 


0.84±0.03 1.08±0.10 

0.80±0.04 0.61±0.04 
0.68±0.04l[5;J]^ 0.52±0.03l[!:2? 


0.46±0.03 0.73±0.06 
0.39±0.02 0.35±0.02 
0.64±0.04 0.46±0.03 
0.54±0.03tg;[5^ 0.32±0.02±0.06 


Hp: EMC @ 
23,4<I^<20 

fiN: E665 Hi 
30,8<1/1^<32 
i^A/^: BBCNC H 
>10,2.5<M^<40 

ep: HI 
300,M^>66 



Table 10: Gaussian fit results for e+e^ collider experiments with > 29 GeV, and from IN 
scattering experiments at high energies. Where two errors are given the first one is statistical, 
the second systematical. All entries refer to fits to the double ratios RR(T) without Coulomb 
and final state corrections. Where available, corrected results are given in brackets. ^ Values 
derived from the single ratios R{T). 
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latter observation confirms, in a new energy regime, earlier work on the search for scale-invariance 
in multi-hadron production. 

Though the evidence for a particle density dependence of the source size, when taken from our 
data alone is marginal, it seems consistent with the trend observed in hadron-hadron collisions. 
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Figure 1: Non-difFractive events - a) Inclusive two-particle densities as a function of T for 
data: full circles denote the like-sign pairs, P2{T); open circles unlike-sign pairs, p2(T); open 
triangles the density obtained from mixing events, pi® pi{T). b) Inclusive two-particle densities 
for reconstructed MEAR Monte Carlo events as in (a). MEAR(BEC) predictions for like-sign 
pairs are indicated as stars, c) Ratios of data to reconstructed MEAR(BEC) for the like and 
unlike-sign densities. 
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Figure 2: a) MEAR Monte Carlo predictions for the invariant mass distributions for unlike- 
sign pion pairs from the decays of the following particles: uj (full circles), p (open circles) and rj 
(shaded-histogram) . The full histogram is the sum of all three contributions. The contribution 
to the like-sign density from rj' decays is shown as stars, b) MEAR Monte Carlo prediction for 
the ratio of densities for vr^vr^ pairs to that of all like-sign charged tracks pairs. 
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Figure 3: Comparison of the non-difFractive data (filled circles) with the MEAR (solid line) and 
MEPS (dashed line) Monte Carlo simulations for the ratios: a) i?""'(r), b) i?'"*(r) and c) 

i?'"(r). 
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Figure 4: Comparison of the difFractive data (filled circles) with the RAPA (solid line) and 
RAPP (dashed line) Monte Carlo simulations for the ratios: a) i?""*(r), b) i?'"*(r) and c) 

i?'"(r). 
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Figure 5: Comparison of the non-diffr active data (filled circles) with the MEAR(BEC) simulation 
for the ratios: a) i?"™(T) and b) B}'"-{T). 
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Figure 6: a) The non-difFractive data for the double ratio RB}^. b) The non-diffractive data 
for the double ratio RR^^. The data points marked by open circles have not been included in 
the fit. The solid curve results from a fit to the Gaussian parametrisation of Eq. 
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Figure 7: The diffractive data for the double ratio (a) and (b). The data points 

marked by open circles have not been included in the fit. The solid curve results from a fit to 
the Gaussian parametrisation of Eq. ([T^). 
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Figure 8: The double ratio for the non-diffractive data sample in intervals of observed 

multiplicity. The curves are fits to the Gaussian parametrisation of Eq. (|14|) . 
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Figure 9: The double ratios a) RR^'^{T), b) RR^''{T) for the non-diffractive data fitted to the 
Gaussian (Eq. (p!^)), exponential (Eq. (p^)) and power law (Eq. ^) parametrisation. 
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Figure 10: The radius parameter r versus charged particle density. The full circles are from this 
analysis, the open circles (triangles) are from pp collisions at a centre of mass energy of 1.8 
TeV |^5| and 630 GeV |^^. Only statistical errors are shown. 
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